CERN-TH/2003-092 
HIP-2003-23/TH 



Q-ball dynamics from atomic Bose-Einstein condensates 

K. Enqvist^ 1 and M. Laine b ' 2,3 

^Dept. of Physical Sciences, University of Helsinki, and Helsinki Institute of Physics, 
P.O.Box 64, FIN-00014 University of Helsinki, Finland 

b Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

Abstract 

Relativistic scalar field theories with a conserved global charge Q possess often (meta)stable 
spherically symmetric soliton solutions, called Q-balls. We elaborate on the perfect formal 
analogy which exists between Q-balls, and spherically symmetric solitons in certain non- 
relativistic atomic Bose-Einstein condensates, for which the dominant interatomic interaction 
can be tuned attractive. In a harmonic trap, present in existing experiments, the Q-ball 
solution is modified in an essential way. If the trap is significantly prolongated in one direction, 
however, then genuine solitons do appear, and actual experimental data can be obtained for 
some of the Q-ball properties studied numerically in the relativistic cosmological context, 
such as their formation and collisions. We also suggest conditions under which the same 
cosmologically relevant analogies could be extended to the fully three-dimensional case. 
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1. Introduction 



Solitons, (meta)stable non-dispersive bound states, appear in all areas of physics described 
by non-linear field theories. They can be considered as lumps of spatially localised matter. 
As is well-known, solitons can be topological or non-topological in nature. Here we shall be 
concerned with non-topological solitons, which are defined as field configurations having the 
same boundary conditions at asymptotic infinity as the vacuum state. For an extensive review 
on non-topological solitons in both relativistic and non-relativistic field theories, see pQ. 

Most soliton solutions exist in one space-dimension only, as required by a theorem due to 
Derrick The theorem can, however, be circumvented, either by including gauge fields of 
non-zero spin or, as will be discussed here, by considering time-dependent solutions. The 
bottom line is that stable spherically symmetric non-topological solitons H], also called 
Q-balls [3], can exist if the theory contains scalar fields with suitable self-interactions, a 
conserved particle number, or "charge" Q, and the charge is carried by massive particles. 

In particle physics and cosmology, Q-balls have attracted a lot of attention recently. The 
reason is that the requirements for their existence are satisfied by (approximately) supersym- 
metric theories fT, considered among the best alternatives for physics Beyond the Standard 
Model. Indeed, such theories do have new scalar fields, in the form e.g. of "squarks". In 
this case the conserved charge is the baryon number B (or some combination of the baryon 
and lepton numbers). The precise properties of Q-balls depend on the particular model of 
supersymmetry breaking, but for many conceivable alternatives, supersymmetry-based Q- 
balls may contribute significantly to the dark matter [7] and baryon contents [SJ EH of 
the Universe (for recent reviews, see Stable Q-balls can also be directly searched for in 

existing and planned experiments 12 . 

In the present paper we transport the Q-ball formalism to non-relativistic atomic Bose- 
Einstein Condensates (BECs), which during the last few years have been the subject of 
exciting experimental developments. Our purpose is to elaborate on the formal analogy that 
exists between Q-balls and various solitons in BECs (see PJEI] for reviews). We thus discuss, 
on one hand, whether it might be possible to observe spherically symmetric three-dimensional 
(3D) Q-balls in actual BEC experiments and, on the other hand, what kind of analogies can be 
drawn from the already existing experiments with essentially one-dimensional (ID) solitons, 
to support numerical studies of their dynamics. Our hope is that these analogies might allow 
to obtain some insights on the behaviour of Q-balls also in cosmology. 

2. Q-balls in relativistic field theory 

In order to set up the relation to the non-relativistic case, let us start by briefly reiterating 
the properties of non-topological solitons, or Q-balls, at zero temperature in a relativistic 
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field theory (for a recent review see, e.g., ^H) 4 . Consider a generic field theory containing 
a scalar field <j) and having a global U(l) symmetry. Let us denote the scalar potential by 
C/(|(/)|). The Minkowskian action is then 



Sm 



dt d 3 x 



(2.1) 



According to the Noether theorem, the system described by this action possesses a conserved 
"charge", Q, 



Q 



d s x 



i(fd t ^-(j)dt4>* 



(2.2) 



In the non-relativistic case Q could be the number of atoms, while in the relativistic case it 
could be the baryon number, as in super symmetric theories. 

The question is, what kind of solutions are there for the classical equations of motions 
derived from Sm, given some fixed value of Q? In order to answer this question it is convenient 
to introduce a Lagrange multiplier \x conjugate to Q, and consider the expression for the 
energy of the system in the sector of a fixed \x first 5 (see, e.g., E3)- Then, a Q-ball 
solution has the form @] 

<fr(x.,t) = exp(— ifj,t)(f)(x.). 



The energy related to this configuration in the given ensemble is 



d 3 x 



The chemical potential /i is related to the total charge of the solution Q through 



Q 



dfj, 



2fi J d 3 x | 



(2.3) 



(2.4) 



(2.5) 



Finally, the energy in the sector of a fixed charge is obtained by a Legendre transform, 

E(Q) = + [iQ , (2.6) 

where fi is expressed in terms of Q by inverting Eq. (|2.5j) . 
To study whether 0(/i) has non-trivial extrema, one writes 



(2.7) 



A spatial variation in a(x) costs energy, so that we may assume it a constant, and without loss 
of generality, choose it to vanish. Moreover, energy is minimised by a spherically symmetric 



4 In this section we employ the conventional natural units h = c = 1. 

5 Oftentimes the notation lu is used instead of but in a thermodynamic sense the quantity in question is 
precisely the chemical potential for Q, which is why we prefer this notation. 
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configuration [T^]. The corresponding profile v(r) is then determined by the classical equation 
of motion following from 

( d 2 2 d \ . , 2 dU(v/V2) , . 

By choosing a suitable fi this system does have a non-trivial localised solution, called a Q-ball, 
provided that U(v/\^2)/v 2 has a minimum for some w > 0HIS]. Physically, this implies 
the existence of an attractive interaction which can bind a collection of elementary quanta 
into a condensate, or "lump", of semiclassical matter. 

The fact that a solution exists for the classical equation of motion does not yet guarantee 
that it is stable with respect to quantum fluctuations. This is the case only if the energy 
of the solution, E(Q), is smaller than the energy of an ensemble of free particles of mass m 
carrying the same charge: 

E(Q) < mQ . (2.9) 

Under these conditions Q-balls are absolutely stable [I]- If Eq. 1)2. 9j) is not satisfied, then 
Q-balls can still be metastable [I] but possibly long-lived (see, e.g., [T7|). 

These basic considerations can be refined in a number of ways. For instance, finite tem- 
perature corrections can be addressed through the grand canonical potential, Q(T, V,fi). It 
is a standard procedure to derive a Euclidean ("imaginary time") path integral expression 
for Q(T,V,n), and one can generically carry out also "dimensional reduction" in this ex- 
pression, integrating out the non-zero Fourier modes for the dependence of the fields on the 
time-like coordinate. The result is just Eq. (|2.4[) . only with modified parameters, containing 
now all relevant dependence on the temperature T jUJ. One can also address a wide variety 
of different potentials: for instance, if the potential is "flat" at large |0|, modulo possible log- 
arithmic corrections, then the energy of the solution E(Q) scales as E(Q) ~ MQ 3 / 4 01 CB], 
allowing to satisfy Eq. (|2,9j) for Q 3> (M/m) 4 , and making Q-balls absolutely stable. With 
other potentials the growth may be slower than |</>| 2 only by radiatively induced logarithmic 
corrections but it is still possible to find regions in the parameter space where Q-balls 
are absolutely stable [20] • Finally, it is possible to address the formation of Q-balls from 
the fragmentation [5J of an essentially homogeneous initial condensate [221 as wen 
as, in case Q-balls are only metastable, their decays and lifetime ^7|> particularly at finite 
temperatures [71 IHIFH ITU] . 



3. Q-balls in non-relativistic field theory 



Let us now turn to the non-relativistic case appropriate for atomic BECs. As is conventional 
in this context, we denote the scalar field by t/j instead of cf>, and reintroduce Ti. 
The "vacuum" action describing the weakly interacting atoms can be written as 



Sm 



dt d 3 x\ ihi/}*d t ip 



L^ 2 



V + V(x) 



2irh 2 a, 



in 



(3.1) 
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Here m is the atom mass, a is the s-wave scattering length, and V(x) is a possible external 
potential. The conserved Noether charge corresponding to this theory is just the particle 
number, 

Q=/d 3 x|^| 2 , (3.2) 



and the energy for stationary configurations, obtained for an ensemble with a chemical po- 
tential [i conjugate to Q, is then 



HO*) = / d 3 x |^|VVf + h + V(x)] H 2 + + ... j . (3.3) 

J { 2m L J 777, J 

As is conventional, the chemical potential has been additively redefined here such that it 
contains the part Vo ~ m in Eq. ()3.1j) . The energy for a fixed charge is again obtained from 

E(Q) = Q(p) + nQ , (3.4) 

where /i is expressed in terms of Q by inverting Eq. (|3.2j) for a given solution ip depending on ji. 
Let us remark that finite temperature effects could be taken into account in complete analogy 
with the relativistic case: one can again write down a Euclidean path integral expression for 
the grand canonical potential and carry out dimensional reduction, to arrive at an expression 
of precisely the form in Eq. (|3.3j) . only with modified parameter values |23j . 

The equation of motion following from n(/x) is the (stationary) Gross-Pitaevskii equa- 
tion ;23|> 

r h 2 _ ATrh 2 n 

ip = . (3.5) 



-— V 2 - n + V x) + |V| 2 + .. 

2m m 



We observe that (apart from trivial changes) Eqs. (|3.3j) . ()3.5|l describe precisely the same 
physics as Eqs. (|2.4|) . (|2.8|) . if V(x) = 0. Note that in the relativistic case, f/(|0|) includes 
also a term quadratic in 4>, U(\(/>\) = m 2 \4)\ 2 + ... , so that m 2 — /j 2 el ~ — 2m fi non _ vc \. 

Since Eqs. (|2.4|) . (|3.3|) are equivalent for the homogeneous case V(x) = 0, the conditions for 
the existence of (meta)stable Q-ball solutions are also equivalent: the potential needs to grow 
more slowly than \ip\ 2 . This can be achieved if there is an attractive interaction between the 
atoms or, equivalently, if the s-wave scattering length is negative, a < 0. This is indeed the 
case for instance for the alkali vapour 7 Li |25| . More generally, the magnitude of a in BECs 
can be tuned in a wide range, including both positive and negative values, using a magnetic 
field close to a so-called Feshbach resonance (see, e.g., EH), as nas been demonstrated also 
for 23 Na [22, 85 Rb |2H], and 133 Cs (22]. In the following, we thus assume that a < 0. 

Obviously, setting just a < in Eq. (|3.3|) is somewhat discomforting, because the theory 
is then not well-defined, being unbounded from below. This implies that the system tends 
to undergo a phase transition to the true ground state, possibly a Bose liquid miED]. It 
is observed experimentally, however, that at least on short enough time scales a weakly 
interacting gaseous phase is still present, even when a < 0. 
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The theory in Eq. ()3.3)) can be explicitly stabilized, however, by adding higher order oper- 
ators, for instance 3 1 

Sn(fj,) = J d 3 x ^A\V{i)*^)\ 2 + S(^» 3 | , (3.6) 

where A parameterises the effective range of the two-body scattering problem, and B the 
amplitude for three-body collisions. Relativistic Q-ball solutions in the case that a is negative 
but B is non- vanishing, have been discussed in [321 • On the other hand, there is a range of 
chemical potentials where we are in the region of the "thick-wall approximation" , and any 
stabilising terms, such as B, can be neglected ^H]- I n the following we will for simplicity 
ignore A but keep B, in order to understand when effects from operators such as those 
in Eq. ()3.6j) are important. Note that in the dilute and (almost) homogeneous limit the 
operator multiplied by B can be argued to be parametrically more important than that 
multiplied by A |3T] . 



4. Solution in homogeneous space 

Let us now consider in more detail the non-relativistic but homogeneous case, that is V(x) = 
in Eqs. (J3.1|) . <|3-3|) . ()3.5|) . but B ^ in Eq. (|3.6f) . The solution resembles very much the 
relativistic one discussed in [32], the main difference being in the relation of ji and Q, but 
for completeness and since the solution does not appear to be widely appreciated in the 
atomic BEC literature, we briefly present some of its main features here, using the notation 
conventional in that context. 



4.1. Equations of motion 

As in Eq. (|2.7j) . we can write the solution of Eq. 1)3.5)) in the form 



(4.1) 



where v > and a can be chosen to vanish. Eq. ()3.5)) then takes the form corresponding 
to Eq. J27 



rr f d 2 2_d 

2m \dr 2 r dr 



2nfi a 9 3 „ , 

H + v 2 + - Bv 4 

m 4 



v = 



(4.2) 



The boundary conditions are that 

v'(0) = 0; lim v(r) = lim v (r) = . (4.3) 

r^oo r— »oo 

Since the system is over-constrained, solutions are only found for specific values of v(0). 
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The parametric dependences of all the properties of the solution of Eqs. (|4.2|) . (|4.3|) can 
easily be found out. A non-trivial solution exists for /U < 0, and in the attractive case we are 
interested in here, a < 0. We can then rescale 



with 



r = fa r , v = va v , B = as , 

'vr?i 2 |a|\ 2 1 



a r = 



2m\n\ 



a„ = 



/ m\fM\ 
2TTh 2 \a\ 



a B = 



rn 



I M I 



(4.4) 



(4.5) 



whereby a common factor can be dropped out from the equation. After this rescaling, 
Eq. (|4.2j) becomes 



d z 2d. , , , 

— + -— )v = v- v 6 + — j3v 



dr 



f dr 



16 



(4.6) 



Given the solution with the boundary conditions corresponding to Eq. (|4.3[) . we can compute 
the dimensionless equivalents of Eqs. (|3.2|) . 1)3. 3 j) . (j3.4j) : 



Q/3 



(3~* 



(32 



oo Y 

dr Airr 2 ■ — v 2 
2 



df4irr 2 ■ 



1 ,rr*l2 1 ~2 

- 2 + - W 2 



1 .4 , 
- v H v 

4 32 



(4.7) 

(4.8) 
(4.9) 



Factors of (3 have been chosen such that rescalings back to physical units contain no |^|'s, 
other than implicitly inside the (3 in Qp,Ep: 



E{Q) 



/mB\i/2, 1 \ 
Us 2 J \2n\a\) 



1 \ 1/2 ft J 



V2mB/ 



4m 



E73 . 



(4.10) 
(4.11) 



4.2. Exact numerical solution 

It is well-known (cf. ref. @]) that a non-trivial solution exists for Eq. (|4.6[) with the boundary 
conditions of Eq. ()4.3|) . for /? < 1. (Because we want the theory to be stable, we also assume 
(3 > 0.) We may recall that the simplest way to understand this is to think of f in Eq. (|4.6[) as 
a time variable, u as a position, and to note that Eq. (|4.6|) then corresponds to the movement 
of a classical particle in a potential P(v) = — \ v 2 + j v 4 — ^v 6 , under the influence of some 
friction. The situation is illustrated in Fig. ^left), for (3 = 0. A solution can be found by a 
simple overshooting-undershooting algorithm, and is also illustrated in Fig. ^ We find the 
range 

-0(0) = 4.3374. ..2.0 , (4.12) 
for (3 = 0.0 ... 1.0, respectively. The corresponding Qp, Ep are shown in Fig. |21 
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10.0 



Figure 1: Left: the Q-ball solution for the normalised (inverted) potential P(v) = —v 2 /2 + 
v A /4, corresponding to Eq. (|4.6j) with f3 = 0. The dots are equidistant in f, corresponding to 
uniform time intervals in the classical mechanics analogue. Right: the corresponding profiles 
v(r), for = 0.0,0.5,0.825. 



4.3. Analytic considerations 

Let us next discuss some analytic estimates in the regimes of small and large /3. The regime 
of small (3, or small B, is called the "thick-wall" regime: there the Q-ball resembles a lump 
without a separate interior and a boundary; in other words, the boundary (or "wall") is as 
thick as the radius. Some properties of the solution, such as its behaviour at large f, have 
previously been discussed, e.g., in pJIHH]. The quantity Qp behaves as Qp = Q/fi 1 / 2 , where 
Q m 9.4486. From Eqs. (|4.4|) . (|4.5j) . (|4.10j) . this gives the relation 

Inserting into Eq. (|4.4[) . the physical distance is obtained from f as 

, Q 

r = 4m-\a\r ■ — , (4-14) 

Q 

and the particle number density at a given distance is given by 

n(r) = M\r) = l -v 2 (r) • (gV f-L-V . (4.15) 



2 \QJ V47r|o| 

It is observed that the central density is smaller for larger particle numbers, but the size of 
the solution is larger. The solution remains weakly interacting even in the center, 

n(r)\a\ 3 <C 1 , (4.16) 
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Figure 2: Left: the rescaled charge, Qp (Eq. (|4.7|) ). as a function of (3 (Eq. (|4.4[0 . The two 
branches of solutions, obtained with small and large [3, are well visible. Right: the rescaled 
energy of the Q-ball, Ep (Eq. (|4.9|0 . as a function of Qp (Eq. t|4. 7j) ^ . The upper branch 
(corresponding to (3 < 0.14) leads to unstable solitons, while the lower branch (/? > 0.14) 
leads to metastable (Ep > 0) or stable (Ep < 0) ones. 



for large values of Q. Note that all dependence on B has cancelled in Eqs. (|4.L'-}j) - (j4.16[) . 

These thick-wall solutions of the equations of motion are not stable, however. Using Q = 
—d£l(ij)/dn and Eq. (|3.4j) . one can derive E(Q) from the Q(\fi\) in Eq. (|4.13l) . to obtain 

Note that in contrast to the relativistic case, this energy does not contain the particle rest 
masses. The fact that the total binding energy E(Q) is positive, implies that this solution is 
an excited state. In fact, it is not even metastable: following [I], one can inspect Q-conserving 
field variations around the solution v(r), and find that there is a direction in the field space 
where even a small variation leads to a decrease of E{Q). Therefore, the thick- wall solutions 
correspond to an unstable branch, as shown in Fig. |^1 There is, however, another solution 
with the same charge but a lower energy, to which we now turn. 

The so-called thin-wall regime is obtained as (3 = |/i|-B(m/(7r?i 2 |a|)) 2 approaches unity. At 
the same time, v(0) approaches 2.0. In this limit the core of the soliton is essentially in 
a homogeneous phase, corresponding to the global minimum of the theory, and has a well- 
defined boundary, or "wall", which is thin compared with the radius. Following this 
solution could also be called a droplet of Bose liquid. In this limit, clearly, physics depends 
in an essential way on the value of the stabilising terms, in our case the coefficient B. 
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The properties of the thin-wall solution can again be found in the standard way. One may 
first compute the interface tension of a planar wall at /3 = 1, a = Jq dvy/2[—P(v)] =2/3. 
Then Qp may be approximated as a sum of a surface term, 4ttR 2 &, and a volume term, and 
extremising this expression allows to solve for the radius R. Consequently, Qp and Ep are 
easily obtained, to leading order in 1 — (3. We find the charge 

-2\4 1 
W (1^/3) 



M3) ■ ^ 



and the radius „ „ 

V 2 / TT%\a\ V 8vr / v ; 

The binding energy becomes 

£(Q) = -(^V§ + 0(Q 2 / 3 ), (4.20) 
\ m J B 

and, being negative for large Q, the solitons are absolutely stable. At the same time, the 
central density becomes 

^|a| 3 ~ TT , (4-21) 

mis 

independent of Q. Thus the interactions are no longer weak, for a small -B. 

To summarise, soliton solutions exist independent of the value of B. They come, however, 
in two branches, and the branch which remains there in the limit B — » (thick- wall, or 
small f3) corresponds to unstable Q-balls. Therefore stabilising terms, for instance of the 
form in Eq. (|3.6|) . are essential for the properties of stable atomic Q-balls, just as they are 
in the relativistic case. The stabilising terms tend to lead, however, to strong interactions in 
the interior of the soliton, which may in fact resemble a Bose liquid rather than a dilute gas. 

Finally, let us mention that according to Fig. |2J the Q-ball solutions have a minimum 
charge, Qp 41, corresponding in physical units to Q m i n ~ 41(mB / (2h 2 )) 1 ^ 2 (2ir\a\)~ 2 . For a 
vanishing stabilising term, therefore, Q m i n — * 0, and the Q-ball could have any charge. Note, 
however, that quantum corrections become important for small Q [31], and our classical 
analysis is no longer trustworthy. On the side of large Q, it has been suggested [3H] that 
there can also exist a maximal charge, Q max , beyond which the system undergoes a phase 
transition to the stable ("Bose liquid") phase. Whether this can happen depends on the 
ensemble: in our case (large volume, fixed Q) it cannot, because the charge density in the 
stable phase is so large (cf. Eq. (|4.21|) ) that there are simply not enough particles present to 
fill the whole volume with this phase: the Q-ball solution is in fact the optimal configuration, 
and absolutely stable. 



5. Solutions in harmonic traps, and experimental data 

In experiments where atomic BECs are studied, space is not homogeneous, but there is a 
harmonic trap, characterised by the potential V(x) in Eq. (|3.5|h This modifies the solution 
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in a qualitative way. We reiterate here the situation for spherically ( "3D" ) and axially ( "ID" ) 
symmetric potentials 6 . We do not present any new solutions but simply point out how the 
Q-ball picture fits the qualitative pattern of the condensate behaviour in these traps. 

3D case. Condensate solutions for a (nearly) spherical trap, V(x) = muj 2 r 2 /2, were ob- 
tained in j^H]. Their essential properties are as follows. The trap has a finite width, charac- 
terised by l r = \Jfij (muj r ), where uo r could typically be in the range of 100 Hz or so. Because 
of a finite l r , the radius of the soliton cannot grow freely, but is restricted. Therefore, as more 
particles condense, the only way to accommodate them is to increase the central density. This 
behaviour is opposite to either of the branches discussed in Sec. 14.31 Therefore the trapped 
"soliton" belongs to yet a different "branch" of solutions than the genuine Q-balls. Another 
way to express the issue is that in the genuine Q-ball solutions \i < (cf. Sec. 0J, while in 
the trap solution fi > |3l?] . 

Because of the growth of the central density, trapped solutions with large charges are 
unstable. Indeed, once the central density increases beyond a certain limit, various losses be- 
come overwhelming [37], and the condensate collapses, as is also observed experimentally |38j . 
After the collapse, the condensate may start to grow again, only to experience yet another 
collapse later on [HIl]- The collapse happens when Q ~ 0(l r /\a\), imposing an upper limit on 
the charge, or particle number, in the condensate. 

ID case. In the ID case the trap has a small finite width only in two directions, but is 
very long in one direction. It turns out that in this case genuine solitons can be observed. 
The chemical potential corresponding to them is negative, as in Sec. This solution, called 
a "bright soliton", was discussed in detail already in [4"0] . 

The essential properties of the ID solitons can easily be deduced from the Q-ball results 
in Sec. El in the thick- wall limit B — > (see also |42l I43| and references therein). Let us now 
denote by l r the transverse width of the trap, and by a r ,a v the scaling factors in Eq. (|4.4|) . 
In the expression for Q, then, the homogeneous 3D relation Q ~ Q afa 2 gets replaced with 
Q ~ Qa r l 2 .a 2 . Therefore, Q oc l/^l 1 / 2 . Consequently, the radius now scales as r oc a r oc 
1/Q, the central density as v 2 oc a 2 oc |/i| oc Q 2 , the grand canonical potential as 
Sl(yu) oc a r l 2 al\n\ oc |/u| 3 / 2 , and the energy as E(Q) = f2(/x) — ^9f2(/x)/5/x oc — \ |/u| 3 / 2 oc — Q 3 . 
Because the binding energy is negative, these Q-balls are absolutely stable compared with 
the gaseous phase, even for B = 0. The central density grows with particle number, however, 
which may still cause an instability related to the practical experimental setup for large 
particle numbers, like for trapped 3D solitons. 

Let us now note that once they have formed |44| 145] . it is possible to study experimentally 

6 We follow standard terminology although it is not without the danger of some confusion: In the "ID" 
case the trap is narrow in two directions, while in the "3D" case it is narrow in three. 

7 In ID solitons exist even for a > 41 . These are sometimes called "dark solitons" , and consist essentially 
of a kink solution, where the field is zero inside an otherwise homogeneous condensate. 
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the collisions of such ID solitons |45| . It is very interesting that the collision results are qual- 
itatively similar to what has been found in numerical simulations of the relativistic case |46| , 
viz., that Q-balls with opposite phases repel each other Such similarities may pro- 

vide an exciting opportunity for studying super symmetry-based post-inflationary cosmology 
in the laboratory. 

6. Conclusions 

We have emphasised in this paper that a rigorous formal analogy exists between the non- 
topological solitons, or Q-balls, of relativistic field theories, and three-dimensional solitons 
that could be found in atomic BECs with a negative s-wave scattering length. 

In order to be stable, the three-dimensional BEC soliton requires an additional stabilising 
term, beyond the usual four-point interaction. The precise form of the stabilising term is not 
essential, however: as an example, we have considered three-body scattering |31| . the strength 
of which we denoted by B. Whether a stabilising interaction of precisely this type could be 
obtained in atomic BECs with a negative scattering length, either directly or effectively as 
a consequence of some coupling of the atoms to external fields, remains at present an open 
issue. If it exists, as could be argued from general principles following or if some other 
stabilising mechanism takes over |131 I3L)| (ultimately even the hard core atomic repulsion 
should be sufficient), then the Q-ball is stable against decay into its quanta, the free atoms. 

Without any stabilising term, the stationary three-dimensional soliton still exists, but it 
has a finite lifetime. The lifetime is currently unknown, but it could be determined by solving 
the time-dependent Gross-Pitaevskii equation around the solution we have presented here. 

These considerations have implications on both contexts in which Q-balls may appear. 
In the BEC case, it is an interesting question whether the spherical Q-ball solutions in a 
homogeneous space, which are quite different from the traditional 3D trapped ones discussed 
in Sec. 03 could also be observed experimentally. This is no doubt a challenging task. In 
principle one could attempt to tune the trap frequency to as small a value as possible while 
a > 0, and then tune a negative, to collapse an almost homogeneous BEC into a genuine 
Q-ball soliton. Alternatively one could start with a significantly prolongated trap holding a 
genuine ID soliton, and then slowly decrease the trap frequency in the transverse directions, 
to try and restore spherical symmetry. 

On the side of cosmology, where most of the interest in Q-balls has been in recent years, 
the current understanding of their dynamics is based on solving classical equations of motion. 
In the actual BEC experiments, of course, the system contains also quantum and thermal 
fluctuations, modifying the dynamics. Thus experimental results from BECs can to some 
extent test which features of the dynamics are robust with respect to these fluctuations. The 
existing ID experiments US are very encouraging in this respect, but it would of course 
be even more remarkable if they could be extended to the (almost) homogeneous 3D case, as 
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outlined above. 
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